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Contest Associations have been found between day-to-day particulate air pollution 
and increased risk of various adverse health outcomes, including cardiopulmonary mor¬ 
tality. However, studies of health effects of long-term particulate air pollution have 
been less conclusive. 

Objective To assess the relationship between long-term exposure to fine particu¬ 
late air pollution and all-cause, lung cancer, and cardiopulmonary mortality. 

Design, Setting, and Participants Vital status and cause of death data were col¬ 
lected by the American Cancer Society as part of the Cancer Prevention 11 study, an on¬ 
going prospective mortality study, which enrolled approximately 1.2 million adults in 1982. 
Participants completed a questionnaire detailing individual risk factor data (age, sex, race, 
weight, height, smoking history, education, marital status, diet alcohol consumption, and 
occupational exposures) The risk factor data for approximately 500000 adults were linked 
with air pollution data for metropolitan areas throughout the United States and com¬ 
bined with vital status and cause of death data through December 31,1993. 

/Main Outcome Measure All-cause, lung cancer, and cardiopulmonary mortality 

Results Pine particulate and sulfur oxide-related pollution were associated with all¬ 
cause, lung cancer, and cardiopulmonary mortality. Each 10-pg/m 5 elevation in fine 
particulate air pollution was associated with approximately a 4%, 6%, and 8% in¬ 
creased risk of all-cause, cardiopulmonary, and lung cancer mortality, respectively. Mea¬ 
sures of coarse partide fraction and total suspended partides were not consistently 
associated with mortality. 

Conclusion Long-term exposure to combustion-related fine particulate air pollu¬ 
tion is an important environmental risk factor for cardiopulmonary and lung cancer 
mortality. 
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B ased on several severe air 
pollution events, 1 ' 5 a temporal 
correlation between extremely 
high concentrations of particu¬ 
late and sulfur oxide air pollution and 
acute increases in mortality was well 
established by the 1970s. Subse¬ 
quently, epidemiological studies pub¬ 
lished between 1989 and 1996 re¬ 
ported health effects at unexpectedly low 
concentrations of particulate air pollu¬ 
tion/ The convergence of data from 
these studies, while controversial, 3 
prompted serious reconsideration of 
standards and health guidelines 6 ' 10 and 
led to a long-term research program de¬ 
signed ro analyze health-related effects 
due to particulate pollution, 11 ' 13 In 1997, 
the Environmental Protection Agency 
adopted new ambient air quality stan¬ 
dards that would impose regulatory lim¬ 
its on fine particles measuring less than 
2.5 pm in diameter These new 

standards were challenged by industry 
groups, blocked by a federal appeals 
court, but ultimately upheld by the US 
Supreme Court, 16 

Although most of the recent epide¬ 
miological research has focused on ef¬ 


fects of short-term exposures, several 
studies suggest that long-term expo¬ 
sure may be more important in terms 
of overall public health. 6 The new stan¬ 
dards for long-term exposure to PMjj 
were originally based primarily an 2 
prospective cohort studies, 15,16 which 
evaluated the effects of long-term pol¬ 
lution exposure on mortality. Both of 
these studies have been subjected to 
much scrutiny, 5 including an exten¬ 
sive independent audit and reanalysrs 
of the original data. 17 The larger of these 


2 studies linked individual risk (actor 
and vital status data with national am¬ 
bient air pollution data. 10 Our analysis 
uses data from the larger study and 
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(1) doubles the follow-up time to more 
than 16 years and triples the number 
of deaths; (2) substantially expands ex¬ 
posure data, including gaseous copol- 
lurant data and new PMj, data, which 
have been collected since the promul¬ 
gation of the new air quality stan¬ 
dards; (3) improves control of occupa¬ 
tional exposures; (4) incorporates 
dietary variables that account for total 
fat consumption, and consumption of 
vegetables, citrus, and high-fiber grains, 
and (5) uses recent advances in statis¬ 
tical modeling, including the incorpo¬ 
ration of random effects and nonpara- 
meiric spatial smoo thing components 
in the Cox proportional hazards model. 

METHODS 
Study Population 

The analysis is based on data collected 
by the American Cancer Society (ACS) 
as part of the Cancer Prevention S tudy 
II (CPS-II), an ongoing prospective 
mortality study of approximately 1.2 
million adults.‘“• 1 * Individual partici¬ 
pants were enrolled by ACS volun¬ 
teers in the fall of 1582. Participants re¬ 
sided in all 50 states, the District of 
Columbia, and Puerto Rico, and were 
generally friends, neighbors, or ac¬ 
quaintances of ACS volunteers. Enroll¬ 
ment was restricted to persons who 
were aged 30 years or older and who 
were members of households with at 
least 1 individual aged 45 years or older, 
Participants completed a confidential 
questionnaire, which included ques¬ 
tions about age, sex, weight, height, 
smoking history, alcohol use, occupa¬ 
tional exposures, diet, education, mari¬ 
tal sutus, and other characteristics. 

Vital status of study participants was 
ascertained by ACS volunteers in Sep¬ 
tember of the following years: 1984, 
1986, and 1988, Reported deaths were 
verified with death certificates. Subse- 
quendy, through December 31, 1998, 
vital status was ascertained through au¬ 
tomated linkage of the CPS-II study 
population with the National Death In¬ 
dex.' 5 Ascertainment of deaths was 
more than 98% complete for the pe¬ 
riod of 1982-1988 and 93% complete 
after 1988. 15 Death certificates or codes 


for cause of death were obtained for 
more than 98% of all known deaths. 
Cause of death was coded according to 
the International Classification of Dis¬ 
eases, Ninth Revision (1CD-9). Al¬ 
though the CPS-II cohort included ap¬ 
proximately 1.2 million participants 
with adequate questionnaire and cause- 
of-death data, our analysis was re¬ 
stricted to those participants who re¬ 
sided in US metropolitan areas with 
available pollution data. The actual size 
of the analytic cohort varied depend¬ 
ing on the number of metropolitan ar¬ 
eas for which pollution data were avail¬ 
able Table 1 provides the number of 
metropolitan areas and participants 
available for each source of pollution 
data. 

Air Pollution Exposure Estimates 

Each participant was assigned a met¬ 
ropolitan area of residence based on ad¬ 
dress at time of enrollment and 3-digit 
ZIP code area 10 Mean (SD) concentra¬ 
tions of air pollution for the metropoli¬ 
tan areas were compiled from various 
primary data sources (Table 1). Many 
of the particulate pollution indices, in¬ 
cluding PM},, were available from data 
from the lnhalable Panicle Monitor¬ 
ing Network for 1979-1983 ahd data 
from the National Aerometric Data¬ 
base for 1980-1981, periods just prior 
to or at the beginning of the follow-up 
period. An additional data source was 
the Environmental Protection Agency 
Aerometric Information Retrieval Sys¬ 
tem (AIRS). The mean concentration 
of each pollutant from all available 
monitoring sites was calculated for each 
metropolitan area during the 1 to 2 
years prior to enrollment. 17 

Additional Information on ambient 
pollution during the follow-up period 
was extracted from the AIRS database 
as quarterly mean values for each rou¬ 
tinely monitored pollutant for 1982 
through 1998. All quarterly averages 
met summary criteria imposed by the 
Environmental Protection Agency and 
were based on observations made on at 
least 50% of the scheduled sampling 
days at each site. The quarterly mean 
values for all stations in each metro¬ 


politan area were calculated across the 
study years using daily avetage values 
for each pollutant except ozone, For 
ozone, daily 1-hour maximums were 
used and were calculated for the full 
year and for the third quarter Only (ie, 
July, August, September). While gas¬ 
eous pollutants generally had re¬ 
corded data throughout the entire fol¬ 
low-up period of interest, the particulate 
matter monitoring protocol changed in 
[he late 1980s from total suspended par¬ 
ticles to particles measuring less chan 
10 pm in diameter (PM,„), resulting in 
the majority of total suspended par¬ 
ticle data being available in the early to 
mid-1980s and PM :d data being mosdy 
available in the early to mid-1990s. 

As a consequence of the new PM,, 
standard, a large number of sites be¬ 
gan collecting PM U data in 1999. Daily 
PM 2 , data were extracted from the AIRS 
database for 1999 and the first 3 quar¬ 
ters of2000. For each site, quarterly av¬ 
erages for each of the 2 years were com¬ 
puted. The 4 quarters were averaged 
when at least 1 of the 2 corresponding 
quarters for each year had at least 50% 
of the sixth-day samples and at least 45 
total sampling days available. Measure¬ 
ments were averaged first by site and 
then by metropolitan area. Although no 
network of PMj,j monitoring existed in 
che United States between the early 
1980s and the late 1990s, the inte¬ 
grated average of PM j3 concentra¬ 
tions during the period was estimated 
by averaging the PM i3 concentration For 
early and later periods. 

Mean sulfate concentrations for 1960- 
1981 were available for many cities 
based on data from the lnhalable Par¬ 
ticle Monitoring Network and the 
National Aerometric Database. Recog¬ 
nizing that sulfate was artifaciually 
overestimated due to glass fiber filters 
used at that time, season and region- 
specific adjustments were made. 17 Since 
few states analyzed particulate samples 
for sulfates after the early 1980s, indi¬ 
vidual states were directly contacted for 
data regarding filter use. Ion chroma¬ 
tography was used 10 analyze PM w fil¬ 
ters and this data could be obtained 
from metropolitan areas across the 
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Uni ted States. Filters were collected for 
a single reference year (19901 in the 
middle of the 1982-199B study pe¬ 
riod. The use of quartz filters virtually 
eliminated the historical overestima- 
lion of sulfate. Mean sulfate concen¬ 
trations for 1990 were estimated us¬ 
ing sulfate from AIR5, data reported 
directly from individual states, and 
analysis of archived filters 

Statistical Analysis 

The basic statistical approach used in this 
analysts is an extension of the standard 
Cox proportional hazards survival 


model, 31 which has been used for risk 
estimates of pollution-related mortal¬ 
ity in previous longitudinal cohort Stud¬ 
ies. ul< The standard Cox model implie- 
itly assumes that observations are 
statistically independent after control¬ 
ling for available risk factors, resulting 
in 2 concerns with regard to risk esti¬ 
mates of pollution-related mortality. 13 
First, if the assumption of statistical in¬ 
dependence is not valid, the uncer¬ 
tainty in the risk estimates of pollution- 
related mortality may be misstated. 
Second, even after controlling for avail¬ 
able risk factors, survival rimes of par¬ 


ticipants living in communities closer to¬ 
gether may be more similar than 
participants living in communities far¬ 
ther apart, which results in spatial au¬ 
tocorrelation. If this spatial autocorre¬ 
lation is due to missing or systematically 

mismeasured risk factors that ate spa¬ 
tially correlated with air pollution, then 
the risk estimates of pollution-related 
mortality may be biased due to inad¬ 
equate control of these factors. There¬ 
fore, in this analysis, the Cox propor¬ 
tional hazards model was extended by 
incorporating a spatial random-effects 
component, which provided accurate es- 


Tlbl* 1. Summary erf Alternative Pollution Indices* 

Pollutant 
(Years of Data 
Collection] 

Units 

Source of Data 

Deta 

Compilation 

Teamt 

No. ol 

Metropolitan 

Areas 

No. of 

Participants, 
in Thousand b 

Mean (SD) 

PMss 

WW 






1979-1983 


IPMN 

H£l 

61 

359 

21.1 (4.61 

1999-2000 


AIRS 

NYU 

116 

£00 

14.01301 

Average 




Si 

319 

17.7(3.71 

PMiq 

M9frrr’ 






1902-1998 


AIRS 

NYU 

102 

415 

28.6 (5.9( 

PM,! 

(jg/rp 






1979-1983 


IPMN 

HE! 

63 

359 

40.3 (7.7) 


ug/m 3 






1979-1983 


IPMN 

HS 

63 

399 

19 2 16-1) 

Total suspended particles 

MS/rrt 1 






1980*1981 


NAD 

HEI 

156 

590 

B8.0(16.7) 

1979*1903 


IPMN 

HEI 

68 

351 

73.7 (14.3) 

1982-1899 


AIRS 

NYU 

150 

573 

56.7 (13.1) 

Sulfate 

pgfrrP 






1980-1961 


IPMN Bnd NAD. 
antfacr adjured 

HEI 

149 

572 

6.5 (2.0) 

1990 


Compilation and analysis 
of pm, , fillers 

NYU 

53 

269 

6.2 (2.0) 

Sulfur dioxide 

peb 

AIRS 





1980 



HB 

118 

S20 

9.7 («3) 

1982-1996 



NYU 

126 

539 

6.7 (3.0) 

Nftrogen dksdde 

ppb 

AIRS 





1080 



HB 

78 

<109 

27.9 (92) 

1982-1990 



NYU 

101 

493 

21.4 (7.1) 

Carbon moncaide 

ppm 

AIRS 





1960 



HEI 

113 

519 

17 (0.7) 

1962-1996 



NYU 

122 

536 

1.1 (0,4) 

Ozone 

ppb 

AIRS 





1980 



HB 

13* 

569 

47.9(11.0) 

1982-1996 



NYU 

lie 

525 

45 5(7.3) 

1982-19981 



NYU 

134 

5S7 

59,7 02.6) 


-PM J6 mtfeatea particles measuring leas man 2.5 um m tfamowr, PM*), particles measuring lags irwn TO pm in ciameter: PM* particlee measuring less than 15 um in demaier. 
PM,,**, caridas maaattfir^ berween 2.5 and 15 uni In diameter, pg/m*, mcfttgrame per cubic meter; ppb, Darts Dor Wlbft; ppm. pans pdf fftMOb! IPMN. hfiaiabte Panide 
Mcrihoring Network; AIRS. Aeromotric information Retnar^a) SyHiem lEmfronmonrat Protwaton Agency]; and NAD, National Aoromewc P er sftaae. 
tHB Indicated data wra compiled by the Hoatth cfloctf institute reanaiyeia team, which wn prosneusfy p\i*ahed. t? NYU mCicauos data «wi compiled at the Na* Y»k Urwerafty 
School of Medicine, Nelson institute of Envlronmontai Mocfrcme (KJ. and S.D.T.). 

T Daily 1 -bout rrvpo'rTwms wero usad- VaIusb ware cetaulalad only for tho third quarter [b. July. Aug list, Septambff]. 
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□mates of ihe uncertainty of effect esti- variables were classified as either indi- In. addition to the baseline analysis, 
males. The model also evaluated spa- caior (ie, yes/no, binary, dummy) van- several additional sets of analysis were 
tial autocorrelation and incorporated a ables or continuous variables. Van- conducted First, to more fully evalu- 
nonparametric spatial smooth compo- ables used to control for tobacco smoke, ate the shape of the concentration- 
nent (to account for unexplained spa- for example, included both indicator response function, a robust locally 
tial structure). A more detailed descrip- and continuous variables. The smok- weighted regression smoother 26 (within 
tion of this modeling approach is ing indicator variables included; cur- the generalized additive model fratne- 
I provided elsewhere. 21 rent cigarette smoker, former ciga- work 27 ) was used to estimate the rela- 

The baseline analysis in this study es re tie smoker, and a pipe or cigar smoker tionship between particulate air pollu* 

I timated adjusted relative risk (RR) ra- only (all vs never smoking) along with tion and mortality in the second stage 
rios for mortality by using a Cox pro- indicator variables for starting smok- of model fitting. Second, the sensitiv- 
portional hazards model with inclusion ing before or after age 1S years. The ity of the fine particle mortality risk es- 
of a metropolitan-based random- continuous smoking variables in- rimates compared with alternative mod- 
effects component. Model fitting in- eluded; current smoker’s years of smok- cling approaches and assumptions was 
volved a 2-stage process. In the first ing. current smoker’s years of smok- evaluated. Standard Cox proportional 
stage, survival data were modeled us- ing squared, current smoker’s cigarettes hazards models were fit to the data fil¬ 
ing the standard Cox proportional haz- per day, current smoker’s cigarettes per eluding particulate air pollution as a 
aids model, including individual level day squared, former smokers years of predictOT of mortality and sequen- 
covariates and indicator variables for smoking, former smoker's years of tially adding (In a controlled forward 
each metropolitan area (without pol- smoking squared, former smoker’s ciga- stepwise process) groups of variables 
lution variables). Output from stage 1 rettes per day, former smoker’s ciga- to control for smoking, education, man- 
provided estimates of the metropolitan- rettes per day squared, and the num- tal status, BMI, alcohol consumption, 

• specific logarithm of the RRs of mor- ber of hours per day exposed to passive occupational exposures, and diet, 
j tality (relative to an arbitrary refer- cigarette smoke. In addition, to evaluate the sensitW- 

ence community), which were adjusted To control for education, 2 indica- ity of the estimated pollution effect 
for individual risk factors. The corre- tot variables, which indicated comple- while more aggressively controlling for 
lation between these values, which was non of high school or education be- spatial differences in mortality, a 2-di- 
induced by using the same reference yond high school, were included, mensional term to account for Spatial 
community, was then removed. 13 In the Marital status variables included indi- trends was added to the modelsandwas 
second stage, the estimates of ad- cator variables for single and other vs estimated using a locally weighted re- 
jUStcd metropolitan-specific health re- married. Both body mass index (BMI) gression smoother. The “span" param- 
sponses were related to fine particu- values and BMI values squared were in- eier, which controls the complexity of 
! late air pollution using a linear random- eluded as continuous variables. Indi- the surface smooth, was set at 3 difler- 
I effects regression model. 2 * The time cator variables for beer, liquor, and wine ent Settings to allow for increasingly ag- 
} variable used in the models was sur- drinkers and nonresponders vs non- gressive fitting of the spatial structure, 
i vival time from the date of enroll- drinkers were included to adjust for al- These included a default span of 50%, 
ment. Survival times of participants who cohol consumption. Occupational ex- the span that resulted in the lowest un¬ 
did not die were censored at the end of posure was controlled for using various explained variance in morality rate be- 
the study period. To control for age, sex, indicator variables; regular occupa- tween metropolitanarcas, andtbespan 
and race, all of the models were strati- tional exposure to asbestos, chemicals/ that resulted in the strongest evidence 
fiedby 1-year age categories, sex, and acids/solvents, coal orstone dusts, coal (highest P value) to suggest no re- 
race (while vs other), which allowed tar/pitch/asphalt, diesel engine ex- sidual spatial structure. The risk esti- 
each category to have its own baseline haust, or formaldehyde, and addi- mates and SEs (and thus the confi- 
hszard. Models were estimated for all- tional indicator variables that indi- dence intervals) were estimated using 
cause mortality and for 3 separate mor- cated 9 different rankings of an generalized additive modeling 27 with 
tality categories; cardiopulmonary occupational dirtiness index that has S-Plusstatisticalsoftwarewhichpro- 
(lCD-9 401-440 and 460-519), lung been developed and described else- vides unbiased effect estimates, butmay 
cancer (/CD-9 162), and all others. where. ,T ‘ y Two diet indices that ac- underestimate SEs if there is signifi- 

Models were estimated separately for counted for fat consumption and con- cant spatial autocorrelation and signifi- 

each of the 3 fine particle variables, sumption of vegetables, citrus, and cant correlations between air pollu- 
PMj j (1979-1983), PMi. } (1999- high-fiber grains were derived based on tion and the smoothed surface of 
2000), and PM 3l5 (average). Indi- information given in the enrollment morality. Therefore, evidence of spa- 
vidual level covariates were included in questionnaire. 18 Quintile indicator vari- dal autocorrelation was carefully evalu¬ 
ate models to adjust for various impor- ables for each of these diet indices were atedand tested using the Bartlett iesl 22 
tant individual risk factors. All of these also included in the models. 18 The correlations of residual mortality 
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with distance between metropolitan ar¬ 
eas were graphically examined. 

Analyses were also conducted of 
effect modification by age, sex, smok¬ 
ing status, occupational exposure, and 
education. Finally, models were fit us¬ 
ing a variety of alternative pollution in¬ 
dices, including gaseous pollutants. 
Specifically, models were estimated 
separately for each of the pollution vari¬ 
ables listed in Table 1, while also in¬ 
cluding all of the other risk factor vari¬ 
ables. 


RESULTS 

Fine particulate air pollution gener¬ 
ally declined ih the United States dur¬ 
ing the follow-up period of this study. 
Figure 1 plots mean PMjj concentra¬ 
tions for 1999-2000 over mean PMj.j 
concentrations for 1979-1933 for the 


Figure 1. Mean Fine Particles Measuring 
Less Than 2.5 pm in Diameter (PMj.s) 



Mean PMis concentrations in micro grams per me- 
it n eubed for 1$9$-2000 aw plotted alcn# with con- 
ccn*jatf£3nsfor1^79*13fl3forthe5l meiropolfifin ar¬ 
eas with paked pollution data. The dotted line is a 
reference 45"-equality tine. 


51 cities in which paired data were 
available. The concentrations of PMjj 
were lower in 1999-2000 than in 1979- 
1983 for most cities, with the largestTe- 
duction observed in the cities with the 
highest concentrations of pollution dur¬ 
ing 1979-1983. Mean PM. j levels in the 
2 periods were highly correlated 
(r=0,78). The rank ordering of cities 

by relative pollution levels remained 
nearly the same. Therefore, the rela¬ 
tive levels of fine particle concentra¬ 
tions were similar whether based on 
measurements at the beginning of the 
study period, shortly following the 

study period, or an average of the 2. 

As reported in Table 2, all 3 indices 
of fine particulate air pollution were as¬ 
sociated with all-cause, cardiopulmo¬ 
nary, and lung cancer mortality, butnot 
mortality from all other causes com¬ 
bined. FIGURE 2 presents the nonpara- 
raerric smoothed exposure response re¬ 
lationships between cause-specific 

mortality and PM U (average). The log 

RRs for all-cause, cardiopulmonary, and 
lung cancer mortality increased across 
the gradient of fine particulate matter. 
Goodness-of-llt tests indicated that the 
associations were not significantly dif¬ 
ferent from linear associations (P> .20). 

The fine particle mortality RR ratios 
from various alternative modeling ap¬ 
proaches and assumptions are pre¬ 
sented in FIGURE 3. After controlling for 
smoking, education, and marital sta¬ 
tus, the controlled forward stepwise in¬ 
clusion of additional co variates had litde 
influence on the estimated associations 
with fine particulate air pollution on car¬ 
diopulmonary and lung cancer mortal¬ 
ity'. As expected, cigarette smoking was 
highly significantly associated with el¬ 


Table 2. Adjusted Mortality Relative Risk (RR) Associated With a 10-pg/i»i s Change in Fine 
Particles Measuring Less Than 2,5 jim in Diameter 


Adjusted RH (95% Cl}" 


Cause of Mortality 

1979-1 983 

1999-2000 

1 

Average 

All-cause 

1.04(1.01-1,09) 

1.06 ( 1 . 02 - 1 . 10 ) 

1.06(1.02-1.11) 

Cardiopulmonary 

1.06(1.02-1.10) 

1.00(102-1.14) 

1.09(1.03-1.16) 

Lung cancer 

1 . 06 ( 1 . 01 - 1 . 16 ) 

1,13(1.04-1.22) 

1.14(1.04-1.231 

Ail other cause 

1.01 (0.97-1.05) 

1.01 (0.07-1.08) 

1.01 (0.95-1-06) 

“Estimated and ac^ustgd based on iho baseline random-effects Cox proportion^ ha»rtf$ model, controlling lor 990 , 
sex, mca. anwklng, education, marfai status, body masc. alcohol ccrisummm occupational exposue. and diet 

. O iruScaies contldanoe Interval. 


evated risk of all-cause, cardiopulmo¬ 
nary, and lung cancer mortality 
(PC.OOl). Estimated RRs for an aver¬ 
age current smoker (men and women 
combined, 22 cigarettes/day for 33.5 
years, with initiation before age 18 years) 
were equal to 2.58, 2.89, and 14.80 for 
all-cause, cardiopulmonary, and lung 
cancer mortality, respectively. Statisti¬ 
cally significant, butsiibstanually smaller 
and less robust associations, were also 
observed for education, marital status, 
BM1, alcohol Consumption, occupa¬ 
tional exposure, and diet variables. Al¬ 
though many of these covanates were 
also statistically associated with mortal¬ 
ity, the risk estimates of pollution- 
related mortality were not highly sen¬ 
sitive to the inclusion of these additional 
covariates. 

Figure 3 also demonstrates that the 
introduction of the random-effects com¬ 
ponent to the model resulted in larger 
SEs of the estimates and, therefore, 
somewhat wider 95% confidence in¬ 
tervals. There was no evidence of sta¬ 
tistically significant spatial autocorre¬ 
lation in the survival data based on the 
Barden test (P>.20) after controlling 
for fine particulate air pollution and the 
various individual risk factors. Further¬ 
more, graphical examination of the cor¬ 
relations of the residual mortality with 
distance between metropolitan areas did 
not reveal significant spatial autocor¬ 
relation (results not shown). Never¬ 
theless, the incorporation of spatial 
smoothing was included to further in¬ 
vestigate the robustness of the esti¬ 
mated particulate pollution effect. Effect 
estimates were not highly sensitive to 
the incorporation of spatial smooth¬ 
ing to account for regional clustering 
or other spatial patterns in the data. 

FIGURE 4 presents fine panicle air 
pollution-related mortality RR ratios af¬ 
ter stratifying by age, sex, education, 
and smoking Status, and adjusting for 
all other risk factors. The differences 
across age and sex strata were not gen¬ 
erally consistent or statistically signifi¬ 
cant. However, a consistent pattern 
emerged from this stratified analysis: the 
association with particulate pollution 
was stronger for both cardiopulmo- 
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nary and lung cancer mortality for par¬ 
ticipants with less education. Also, for 
both cardiopulmonary and lung can¬ 
cer mortality, the RR estimates were 
higher for nonsmokers 

Figure 5 summarizes the associa¬ 
tions between mortality risk and air pol¬ 
lutant concentrations listed in Table 1, 
Statistically significant and relatively 
consistent mortality associations ex¬ 
isted for all measures of fine particu¬ 
late exposure, including PM 2 s and sul¬ 
fate particles. Weaker less consistent 
mortality associations were observed 
wich PMj d and PM !} . Measures of the 
coarse particle fraction (PM w . li3 ) and 
total suspended particles were no i con¬ 
sistently associated with mortality. Of 
the gaseous pollutants, only sulfur di¬ 
oxide was associated with elevated mor¬ 
tality risk. Interestingly, measures of 
PM j j were associated with all-cause car¬ 
diopulmonary, and lung cancer mor¬ 
tality, but not with all other mortality. 
However, sulfur oxide pollution (as 
measured by sulfate particles and/or sul¬ 
fur dioxide) was significantly associ¬ 
ated with mortality from all other causes 
in addition to all-cause, cardiopulmo¬ 
nary, and lung canter mortality. 

COMMENT 

This study demonstrated associations be¬ 
tween ambient fine paniculate air pol¬ 
lution and elevated risks of both cardio¬ 
pulmonary' and lung cancer mortality. 
Each 10-pg/m 5 elevation in long-term av¬ 
erage PMjj ambient concentrations was 
associated with approximately a 4%, 6%, 
and 8% increased risk of all-cause, car¬ 
diopulmonary, and lung cancer mortal¬ 
ity, respectively, although the magni¬ 
tude of the effect somewhat depended 
on the time frame of pollution monitor¬ 
ing. In addition, this analysis addresses 
many of the important questions con¬ 
cerning the earlier, more limited analy¬ 
sis of the large CPS-II cohort, includ¬ 
ing the following issues. 

First, does the apparent association 
between pollution and mortality per¬ 
sist with longer follow-up and as the co¬ 
hort ages and dies? The present analy¬ 
sis more than doubled the follow-up 
time to more than 16 years, resulting 


in approximately triple the number of 
deaths, yet the associations between 
pollution and mortality persisted. 

Second, can the association between 

fine particulate air pollution and in¬ 
creased cardiopulmonary and lung can¬ 
cer mortality be due to inadequate con¬ 
trol of important individual risk factors? 
After aggressively controlling for smok¬ 
ing, the estimated fine particulate pol¬ 
lution effect on mortality was remark¬ 
ably robust. When the analysis was 
stratified, by smoking status, the esti¬ 
mated pollution effect on both cardio¬ 
pulmonary- and lung cancer mortality 
was strongest for never smokers vs 
former or current smokers. This analy¬ 
sis also controlled for education, mari¬ 
tal status, BM1, and alcohol consump¬ 
tion. This analysis used improved 
variables to control for occupational ex¬ 
posures and incorporated diet variables 
that accounted for total fat consump¬ 
tion, as well as for consumption of veg¬ 
etables, citrus, and high-fiber grains. The 
mortality associations wich fine particu¬ 
late air pollution were largely unaf¬ 
fected by the inclusion of these indi¬ 


vidual risk factors in the models. The data 
on smoking and other individual risk fac¬ 
tors, however, were obtained directly by 
questionnaire at time of enrollment and 
do notreflect changes that may have oc¬ 
curred following enrollment. The lack of 
risk factor follow-up data results in some 
rmsdassification of exposure, reduces the 
precision of control for risk factors, and 
constrains our ability to differentiate time 
dependency. 

Third, are the associations between 
fine particulate air pollution and mor¬ 
tality due to regional or other spatial dif¬ 
ferences that are not adequately con¬ 
trolled for in the analysis? If there ate 
unmeasured or inadequately modeled 
risk factors lhat are different across lo¬ 
cations, then spatial clustering will oc¬ 
cur. If this clustering is independent or 
random across metropolitan areas, then 
the spatial clustering can be modeled 
by adding a random-effects compo¬ 
nent to the Cox proportional hazards 
model as was done in our analysis. The 
clustering may not be independent or 
random across metropolitan areas due 
to inadequately measured or modeled 


Figure 2. Nonparametric Smoothed Exposure Response Relationship 
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Figure 3. Mortality Relative Risk (RR) Ratio Associated With 10-yg/m 5 Differences of PM J3 Concentration; 
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risk factors (either individual or eco¬ 
logical). If these inadequately mea¬ 
sured or modeled risk factors are abo 
spatially correlated with air pollution, 
then biased pollution effects estimates 
may occur due to confounding. How¬ 
ever, in this analysis, significant spa¬ 
tial autocorrelation was not observed 
after controlling for fine particulate air 
pollution and the various individual risk 
factors. Furthermore, to minimize any 
potential confounding bias, sensitiv¬ 
ity analyses, which directly modeled 
spatial trends using nonparametric 
smoothing techniques, were con¬ 
ducted. A contribution of this analysis 
is that it included the incorporation of 
both random effects and nonparamet¬ 
ric spatial smoo thing components to the 
Cox proportional hazards model- Even 
after accounting for random effects 
across metropolitan areas and aggres¬ 
sively modeling a spatial structure that 
accounts for regional differences, the as¬ 
sociation between fine particulate air 
pollution and cardiopulmonary and 
lung cancer mortality persists. 

Fourth, is mortality associated pri¬ 
marily with fine particulate air pollu¬ 
tion or is mortality also associated with 
ocher measures of particulate air pol¬ 
lution, such as PMib, total suspended 
particles, or with various gaseous pol¬ 
lutants? Elevated mortality risks were 
associated primarily with measures of 
fine particulate and sulfur oxide pol¬ 
lution. Coarse particles and gaseous pol¬ 
lutants, except For sulfur dioxide, were 
generally not significantly associated 
with elevated mortality risk. 

Fifth, what is the shape of the con¬ 
centration-response function? Within the 
range of pollution observed in this analy¬ 
sis, the concentration-response func¬ 
tion appears to be monotonic and nearly 
linear. However, this does not preclude 
a leveling off (or even steepening) at 
much higher leveb of air pollution. 

Sixth, how large is the estimated mor¬ 
tality effect of exposure to fine particu¬ 
late air pollution relative to other risk fac¬ 
tors? A detailed description and 
interpretation of the many individual risk 
factors that are controlled for in the 
analysis goes well beyond the scope of 


Figure 4. Adjusted Mortality Relative Risk (RR) Ratio Associated With 10-(jg/m' Differences 
of PMjj Concentrations 
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this report- However, the mortality risk 
associated with cigarette smoking has 
been well documented using the CP5-11 
cohort. 11 The risk imposed by exposure 
to fine paniculate air pollution is obvi¬ 
ously much smaller than the risk of ciga¬ 
rette smoking. Another risk factor that 
has been well documented using the 
CPS-II cohort data is body mass as mea¬ 
sured by BMI. 50 The Word Health Or¬ 
ganization has categorized BM1 values 
between 16 5-24.9 kg/rrr as normal; 25- 

29.9 kg/m J , grade I overweight, 30- 

39.9 kg/m 2 , grade 2 overweight; and 40 
kg/m 1 or higher, grade 3 overweight. 11 
In the presenL analysis, BM1 values and 
EMI values squared were included in the 
proportional hazards models. Consis¬ 
tent with previous ACS analysis, 30 BMI 
was significantly associated with mor¬ 
tality, optimal BMI was between ap¬ 
proximately 23.5 and 24.9 kg/m 1 , and 
the RR of mortality for diflerentBMI val¬ 
ues relative to the optimal were depen¬ 
dent on sex and smoking s ta tus. For ex- 

i ample, the RRs associated with BMI 
values between 30.0 and 31.9 kg/m 3 (vs 
Optimal) would be up to approxi¬ 


mately 1.33 for never smokers. Based on 
these calculations, mortality risks asso¬ 
ciated with fine particulate air pollu¬ 
tion at levels found in more polluted US 
metropolitan areas are less than those as¬ 
sociated with substantial obesity (grade 
3 overweight), but comparable with the 
estimated effect of being moderately 
overweight (grade 1 to 2). 

In conclusion, the findings of this 
study provide the strongest evidence to 
date that long-term exposure to fine par¬ 
ticulate aii pollution common to many 
metropolitan areas is an important risk 
factor for cardiopulmonary mortality In 
addition, the large cohort and extended 
follow-up have provided an unprec¬ 
edented opportunity to evaluate asso¬ 
ciations between air pollution and lung 
cancer mortality. Elevated fine particu¬ 
late air pollution exposures were asso¬ 
ciated with Significant increases in lung 
cancer mortality. Although potential ef¬ 
fects of other unaccounted for factors 
cannot be excluded with certainty, the 
associations between fine particulate ait 
pollution and lung cancer mortality, as 
welt as cardiopulmonary mortality, are 
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